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Abstract. We discuss recent developments and present new findings in the colloidal 
description of soft polymeric macromolecular aggregates. For various macromolecular 
architectures, such as linear chains, star polymers, dendrimers and polyelectrolyte 
stars, the effective interactions between suitably chosen coordinates are shown to 
be ultrasoft, i.e., they either remain finite or diverge very slowly at zero separation. 
As a consequence, the fluid phases have unusual characteristics, including anomalous 
pair correlations and mean- field like thermodynamic behaviour. The solid phases can 
exhibit exotic, strongly anisotropic as well as open crystal structures. For example, the 
diamond and the Al5-phase are shown to be stable at sufficiently high concentrations. 
Reentrant melting and clustering transitions are additional features displayed by such 
systems, resulting in phase diagrams with a very rich topology. We emphasise that 
many of these effects are fundamentally different from the usual archetypal hard sphere 
paradigm. Instead, we propose that these fluids fall into the class of mean-field fluids. 



PACS numbers: 82.70.Dd, 82.70.-y, 6I.25.Hq, 61.20. Gy 
1. Introduction 

One major advantage of soft matter systems in comparison to their atomic counterparts 
is the abihty one has to engineer the constituent particles at the molecular level. In this 
way, an enormous variety in architectures can be achieved, leading to a corresponding 
richness in the structural and phase behaviour of such systems. Polymers play a 
prominent example within this class of materials. They come in a variety of forms, 
such as, e.g., as linear chains, branched, star-shaped, dendritic, copolymer, as well 
as functions, e.g., steric stabilisers, additives, depletants etc. Additional flexibility 
arises from the possibility of influencing the structural and phase behaviour of polymer 
solutions by changing the solvent quality. 

From the theoretical point of view, the task of bridging the gap between the 
microscopic and macroscopic length scales of polymeric systems is a formidable one. 
Indeed, the constituent macromolecules may contain thousands or even millions of 
atoms, interconnected to one another in complicated ways. The program of starting with 
the individual interactions between monomers and proceeding to the calculation of the 
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free energy of the system rapidly becomes intractable. However, considerable progress 
can be made if one invokes a "coarse-graining" procedure, a point of view that has been 
proved very fruitful in many areas of research in condensed matter physics. Here we 
follow a two-step procedure. First, instead of attempting to carry out the statistical 
trace over all the individual monomers in one step, a certain generalised coordinate 
of the macromolecule as an effective point particle is invoked. Examples include the 
centre of mass of the polymer or some suitably selected monomer, as we explain below. 
All monomers belonging to the macromolecules are then traced out for a given, fixed 
configuration of the effective coordinates, which defines an effective interaction between 
these coordinates |^. Once this is achieved, the second step consists of viewing the the 
macroscopic system as collection of point particles interacting by means of the effective 
interaction. Now all known tools from the theory of atomic and molecular fiuids can 
be employed to derive structural and thermodynamic quantities for the system under 
consideration. 

In the recent years, the program sketched above has been carried out with success 
for various polymeric systems. It has been found that the effective interactions obtained 
belong to a new class of ultrasoft potentials which have very unusual properties when 
compared with the hard-sphere (HS) system, the prototype of atomic liquids. In this 
work, we first present a concise review of these novel properties for the fiuid phases 
and then some new results regarding the rich variety of crystal structures that can be 
stabilised by ultrasoft potentials. The rest of the paper is organised as follows. In 
section |^ we give the general definition of the effective interaction and present specific 
examples for a number of systems that have been worked out recently. In section ^ we 
show that the systems described by ultrasoft interactions are well described by a simple 
mean-field picture in the fiuid phase for a wide range of thermodynamic conditions. To 
further delineate these systems we define two categories of mean-field fiuids, the strong 
and the weak ones. The former are characterised by a direct correlation function c(r) 
that satisfies to excellent accuracy the condition c(r) = — /9f (r) in broad regions of 
the thermodynamics space, where v{r) is the interparticle pair potential. The latter 
only satisfy an approximate mean-field picture for their thermodynamic properties, and 
not for their structure. In section ^ we demonstrate that strong mean-field fiuids can 
be further divided into two categories: those displaying reentrant melting and those 
displaying a cascade of clustering transitions, the criterion being set by the positivity of 
the Fourier transform of the effective interaction. In section ^ we discuss the richness 
of the phase diagrams of weak mean-field fiuids, taking the case of star polymers as a 
concrete example. Reentrant melting as well as a wealth of structural phase transitions 
and exotic crystal phases are all shown to be stemming from the ultrasoftness of the 
effective interactions. Finally, in section ^ we summarise and conclude. 
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2. Effective interactions between polymeric macromolecules 

The effective interaction between fiexible, fiuctuating aggregates can be given a precise 
statistical-meclianical definition. Let us consider a solution containing M polymeric 
macromolecules, each one of them composed of monomers. The total number of 
particles in the system is Ai = M x N. One starts from the full Hamiltonian 7i of the 
problem, assumed to be known. Then, out of the Ai particles in the problem (in our case 
all monomers), one selects the M ones that are to be considered as "effective particles" 
and holds them fixed in some prescribed configuration {Ri, R2, . . . , Rm}, where Rj is 
the position of the i-th effective particle. Afterwards, the A4 — M remaining particles 
are canonically traced out and the result of this integration is a constrained partition 
function (5(Ri, R25 • • • , R-m)- The effective Hamiltonian TCes is defined as: 



where f5 = [ksT)'^, with the absolute temperature T and Boltzmann's constant ks- 
It can be shown [0] that such an effective Hamiltonian has two useful properties: 
it preserves the overall thermodynamics of the system and it guarantees that the 
correlation functions of any order between any of the M remaining particles remain 
invariant, regardless of whether the expectation values are calculated with the original 
Hamiltonian Ti. or with the effective one Ties- Though the procedure of tracing out the 
Ai — M degrees of freedom necessarily generates interactions of all orders between the M 
particles Q in many cases it is sufficient to truncate those at the pair level, introducing 
thereby the pair-potential approximation. The great advantage of employing this point 
of view is that, in comparison with the original problem, the numbers of particles, A4, 
has been reduced by a factor N. In addition, whereas in the original problem the pair 
interactions between the monomers are quite complicated, due to the need of taking into 
account the connectivity and architecture of the molecule, in the effective description 
the pair potential is spherically symmetric, depending only on the magnitude of the 
separation vector between the two effective coordinates. A new picture of the polymers 
emerges thereby in which the latter can be seen as ultrasoft colloids having a dimension 
of the order of their radius of gyration Rg. This sets at the same time the characteristic 
length scale of the effective interaction between them. When the polymers are neutral, 
Rg is the only length scale appearing in this colloidal description, whereas if they carry 
charge, additional scales set, e.g., by the concentration of the solution, the counterions 
and/or the salt ions also come into play. We examine some characteristic cases below. 

Polymer chains. Two possible choices for the effective coordinates have been 
investigated thus far. One possibility is to consider the effective interaction fcom(^) 
between the centres of mass of the linear chains, when these are kept fixed at a distance r 
from one another. This was first done in the pioneering work of Flory and Krigbaum ||^, 
who found a Gaussian interaction between the centers of mass. Although the functional 
form of the Flory-Krigbaum potential is correct, the dependence of the pref actor on the 
degree of polymerisation A^ is not: whereas the Flory-Krigbaum mean- field approach 
predicts a A^^/^-dependence the prefactor turns out to be A^-independent for sufficiently 




(1) 
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large A^, as can be easily seen from examining the polymer second virial coefficient. 
Standard scaling arguments show that the radius of gyration Rg is the only relevant 
length scale for the dilute and semi-dilute regimes of polymers in a good solvent . This 
immediately implies that the second- virial coefficient should scale as B2 ^ Rg, which, 
in turn, implies that effective interaction must have an amplitude that is independent 
of Rg 10], at least in the scaling limit. 



A number of simulational fg, |y, [10| as well as theoretical approaches [[l^] involving 



two self-avoiding chains, reached the conclusion that the aforementioned interaction has 
a Gaussian form. The lack of divergence of this effective interaction at zero separation 
should not be surprising. Indeed, the centres of mass of two polymer chains can coincide 
without any of the monomers violating the excluded volume conditions. In addition, it 
can be seen that the effective "particles" one chooses for the coarse-grained description 
of the system do not need to be real particles of the physical system. 



Recently, Louis et al. [|T2|, |T^, |TJ] have independently carried out state-of-the art 
simulations involving not just two but Nc chains and varying the number of chains to 
cover a very broad range of concentrations, ranging from dilute solutions up to nine 
times the overlap concentration. They confirmed that the effective potential has a 
Gaussian-like form which at zero density can be well approximated by 

v^om{r) = eexp[-{r/a)% (2) 

where e = 1.87 ksT and a = 1.08 Rg. For higher densities a superposition of three 
Gaussians provides a very accurate fit |T^, but the basic shape does not deviate much 
from the low-density Gaussian form. The rather weak density dependence can be shown 



to arise from the density independent many-body forces |16[. Moreover, the same 



authors have shown that employing this effective interaction leads to a very accurate 
description of the thermodynamics (equation of state) of polymer solutions for a wide 
range of concentrations, thus confirming the validity of the idea that polymer chains 



can be viewed as soft colloids [14 



An alternative is to consider the end-monomers or the central monomers of the 
chains as effective coordinates 0. General scaling arguments establish that in this case 
the effective interaction diverges logarithmically with the monomer- monomer separation 
r When the central monomers are chosen, linear chains are equivalent to star 

polymers with / = 2 arms. Motivated by this analogy, Jusufi et al. [1^ derived the 
effective interaction fcm('") by combining monomer- resolved, off-lattice simulations with 
theoretical arguments. The sought-for interaction features in this case a logarithmic 
divergence for small separations, in full agreement with the scaling arguments, and 
crosses over to a Gaussian decay for larger ones: 

r I — forr<as; 

-ksTf^'l , ^^^^ 2rV,2 



18 



^ exp [-r^(r^ - cig)] for r > cTs, 



2rVg 
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where cTs = 0.66Rg and r(/) is a free parameter of the order of l/Rg and is obtained 
by fitting to computer simulation results. For / = 2 the value rcXs = 1.03 has been 
obtained |T^, which, together with the potential in Eq. (j^) above yields for the second 
virial coefficient of polymer solutions the value B2/R^ = 5.59, in agreement with the 
estimate 5.5 < i?2/-Rg < 5.9 from renormalisation group and simulations |n . 

Dendrimers. By employing a simple, mean-field theory based on the measured 
monomer density profiles of fourth-generation dendrimers, a Gaussian function of the 
form (Q) has been shown to accurately describe the effective interaction between the 
centres of mass of these dendrimers |]19|, The prefactor e has in this case a higher 
value then for linear polymers, e = 10 ksT. Small-angle neutron scattering (SANS) 
profiles from concentrated dendrimer solutions are reproduced very well theoretically, 
at least below the overlap concentration c*. 

Star polymers. By chemically anchoring / linear chains on a common core, star 
polymers with functionality / are constructed. In the theoretical analysis of the 
conformations and the effective interactions of stars, the finite size of the core particle 
is ignored, an excellent approximation when the chains are long. The natural choice for 
the effective coordinates is now the position of the central particle, i.e., of the star centre. 
For small functionalities, / < 10, Jusufi et al. |1^ have shown that a logarithmic- Gauss 
potential of Eq. accurately describes the effective interaction. The decay parameter 
T of the Gaussian is /-dependent, for details see Ref. [0. For larger functionalities, 
/ > 10, the Daoud-Cotton |21| blob picture of the stars is valid and the star-star 
interaction potential fss('") reads as p2| , p3[ : 

1 



v^Ar 



-In 



18 ^ 



1 



(Js/r 



[ 1 + 77/2 



exp 



2a. 



for r < (Ts 



for r > (Ts 



(4) 



with the "corona diameter" a. = 0.66 Rg. Both star-star potentials, the one valid for 
/ < 10, Eq. (H) and the one valid for / > 10, Eq. (|), show an ultrasoft logarithmic 
divergence as r 0. The strength of the divergence is controlled by the functionality 
/, so that at the formal limit / ^ oo the interaction (^) tends to the HS-potential. 

Poly electrolyte stars. If the polymer chains of a star polymer contain ionisable 
groups, the latter dissociate upon solution in a polar (aqueous) solvent, leaving behind 
charged monomers and resulting in a solution consisting of charged star polymers and 
counterions. The resulting macromolecules are called polyelectrolyte (PE) stars. In PE 
stars the chains are stretched due to the Coulomb repulsion of the charged monomers. 
The degree of stretching and condensation of counterions on the rods depends on the 
amount of charge and on the Bjerrum length. For moderate to high charging fractions, 
the effective interaction between the centres of the PE-stars have been analysed recently 
by means of computer simulations and theory ||2^, ^ 



This interaction, 



IS 



dominated by the entropic contribution of the counterions that remain trapped within 
the star corona. For a broad range of functionalities and charge fractions, it can be 
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accurately described by the fit: 
CfN, 



c 



1 - 

3(1 



+ 



- 1 



3(1-0 
5(1 + Ka) 



5(1 



— ) expf— K(r 
r / 



^)] 



for r < cr; 



for r > cr, 



(5) 



wliere is tlie number of counterions, a the corona diameter of the PE-star, and k 
the inverse Debye screening length due to the free counterions. Finally, C and C are 
fit parameters, where < C < 1- The last condition ensures that the potential of Eq. 
(§) above tends to a finite value as r — 0. Hence, once more we are dealing with an 
ultrasoft interaction that varies slowly as the particle centres approach one another. 

Ultrasoft interactions therefore describe quite a number of different systems. 
Their common characteristic is that the constituent particles are polymers of various 
architectures that dominate the spatial extent of the aggregates. In other words, one 
expects similar ultrasoft interactions to show up also when one deals with core-shell 
particles, consisting of a solid core and a polymeric shell, whenever the thickness of the 
latter greatly exceeds the radius of the former. In addition, the ultrasoft interactions 
can be tuned by controlling the number of arms, the charge, the length of the chains, 
the generation number (in the case of dendrimers) etc. Hence, it is useful to explore 
the general characteristics of this family of potentials and the ramifications they have 
on the structural and thermodynamic properties of the fluid and crystal phases of such 
systems. 



3. Mean-field fluids 

Motivated by the fact that effective interactions between polymeric colloids can be 
bounded (i.e., finite at all separations r), we examine here in general the properties of 
systems characterised by pair potentials of the form 

f (r) = e(j){r/cr), (6) 

with 0(x) < oo for all x. In Eq. (|^) above, e is an energy scale and a a length 
scale. Moreover, v{r) is non-attractive, i.e., d(f){x)/dx < everywhere. We introduce 
dimensionless measures of temperature and density as 



TC 



TC 



(7) 
(8) 



where fcs is Boltzmann's constant and p = N/V is the density of a system of particles 
in the volume V. We will refer to rj as the 'packing fraction' of the system. 

The key idea for examining the high temperature and/or high density limit of such 
model systems in three and higher dimensions is the following. We consider in general 
a spatially modulated density profile p(r) which does not vary too rapidly on the scale 
a set by the interaction. At high densities, pa^ ^ 1, the average interparticle distance 
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a = p~^/^ becomes vanishingly small, and it holds a <^ a, i.e., the potential is extremely 
long-range. Every particle is simultaneously interacting with an enormous number of 
neighboring molecules and in the absence of short-range excluded volume interactions 



the excess free energy of the system can be cast in the mean-field approximation 
(MFA) to be equal to the internal energy of the system [p^]: 

FM -\j j d'rdVt;(|r - r'|)p(r)p(r'), (9) 

with the approximation becoming more accurate with increasing density. Then, Eq. (|^) 
immediately implies that in this limit the direct correlation function c(r) of the system, 
defined as 

c(|r-r1;p) = - li,„ (10) 

becomes independent of the density and is simply proportional to the interaction, namely 
c(r) = -(3v{r). (11) 



Using the last equation, together with the Ornstein-Zernike relation [2S], we readily 



obtain an analytic expression for the structure factor S{k) of the system as 

S{k) = ^-^ , (12) 

I + pt-^(t){ka) 

where 0(g) = J d^xexp(— iq ■ x)0(x) is the Fourier transform of 0(x). 

Bounded and positive-definite interactions have been studied in the late 1970s by 



Grewe and Klein p9| , |30| . The authors considered a Kac potential of the form: 

i;(r)=7X7r), (13) 

where d is the dimension of the space and 7 > is a parameter controlling the range 
and strength of the potential. Moreover, tpi^x) is a nonnegative, bounded and integrable 
function. Grewe and Klein showed rigorously that at the limit 7 — > 0, the direct 
correlation function of a system interacting by means of the potential ( |T3|) is given by 
Eq. (0) above. The connection with the case we are discussing here is straightforward: 
as there are no hard cores in the system or a lattice constant to impose a length scale, 
the only relevant length is set by the density and is equal to p'"^/^ in our model and by 
the parameter 7"^ in model (|l^). In this respect, the limit 7 — in the Kac model of 
Grewe and Klein is equivalent to the limit p — *• 00 considered here. 

Although the limit of Grewe and Klein corresponds to t — 00 and p — > cxd, 
the relation (0) has been shown to be an excellent approximation at arbitrarily low 



temperatures for high enough densities [ETI and for temperatures t ^ 1 practically at 



all densities |^T|. Hence, the mean-field approximation is valid in a vast range of the 
thermodynamic space of such systems, which has led to their characterisation as mean- 
field fluids (MFF) |]32[. Associated with the structural relation (0) are scaling relations 



of thermodynamic quantities, arising from the compressibility sum rule [2S]: 



/"(p) = -c{k = 0; p) = -An / r'^c{r] p) dr, (14) 
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where /(p) = [3F^^{p) /V , and the primes denote the second derivative. From Eqs. (pi]) 
and ( |I^ ) it then follows 

with ?}(0) = A-n ' r'^v{r) dr. This simple scaling is not at all equivalent to a second 
virial theory. In fact, simple virial expansions have a rather small radius of convergence 
for mean field fluids |I3| , |32| . It then follows that the excess pressure, chemical potential 
and compressibility satisfy the scaling relations Pex ~ P^, Pex ~ P and Xex ~ 
p!3| , p7|.|| All these stem from the validity of the strong structural-MFA relation (pip 
which guarantees the validity of the weaker thermodynamicalMFA relation (0). In what 
follows, we will argue that many ultrasoft systems can still satisfy the thermodynamic 
relation ( |l5l) approximately, without satisfying Eq. (0). To distinguish between the two 
classes, we now qualify the term "mean-field fluids" and call the ones for which both 
the structural and thermodynamic MFA work well strong mean-field fluids. For such 
systems, for which the potentials are typically also bounded, the density functional of 
Eq. (H) has been extended to mixtures [0 allowing a straightforward and transparent 
analysis of the phase separation in the bulk, interfacial |Q and wetting properties 
of such mixtures. 

Let us now then turn our attention to interaction potentials such as those given 
in Eqs. (§) and (^. These diverge at the origin slowly enough, so that the three- 
dimensional integral 

d^rv{r)=4:TT / r^f (r) dr =v{0), (16) 



is finite and equal (by definition) to the value of the Fourier transform of the potential 
v{k) at = 0. It is now impossible to satisfy the strong mean-field condition of Eq. (PH) 
everywhere. Indeed, the direct correlation function c(r) has to remain finite at r = 0, 
whereas the pair potential diverges. Hence, as shown in Fig. |l|(a), there will always exist 
a region in the neighborhood of the origin in which Eq. (0) is violated. At the same 
time, it can be seen in this figure that the extent of this region shrinks with increasing 
density, hence the fluid becomes more 'strong mean-fleld'-like as it gets denser. The 
smaller v{0), the lower the density at which the MFA for S{k), Eq. ([I^), becomes a 
reasonable approximation. 

The discrepancies between c(r) and —(3v{r) become much less important when we 
turn our attention to the thermodynamics. To obtain the excess Helmholtz free energy, 
one needs only the integral of r^c(r), see Eq. (0). As demonstrated in Fig. |l|(b), upon 
multiplication with the geometrical factor r^, the deviations of c(r) from —(3v{r) become 

X Note that for polymers in a good solvent Pex ''^ p^f/CSiz-i) ^ ^2.3 -j-j^g semi-dilute regime. The 
density dependence of the pair-potentials is necessary to properly describe this correction to simple 
MFF behaviour jlj]. But this in turn implies that the extra factor 0.3 in the scaling of the pressure 
arises from the many-body interactions, since these are what cause the density-dependence in the first 
place 111. 



Exotic fluids and crystals of soft polymeric colloids 



9 




Figure 1. (a) Comparison between the direct correlation function of a / = 32 star 
polymer fluid at various densities, obtained by solving the Rogers- Young closure, with 
the mean-field result, —(3v{r). (b) Same as in (a) but for the quantities r^c(r) and 
—(3r'^v{r). 



suppressed, so that we can write, to a very good approximation: 



drr c{r; p) — — drr (3v{r). (17) 
Jo 

Eq. ( [T7| ) together with Eq. (|Hp yield an approximate scahng of the excess free energy of 
the weak mean-field fluids with density that is identical to that of the strong mean-field 
fluids, Eq. (|15D . The accuracy of the approximation for the star polymer fluid with 
/ = 32 arms [Eq. (|)] is shown in Fig. |^. The line labeled as exact free energy there 
was obtained by solving the Rogers- Young closure for the fluid at a wide density 
range and subsequently utilizing the compressibility sum rule [Eq. (|l5)] to obtain the 
excess free energy. Comparisons with simulations have indeed demonstrated that 



this procedure delivers an essentially exact numerical result. 

Clearly, the mean-field approximation improves with increasing density, as the 
number of particles effectively interacting with one another grows. The crossover density 
Px above which the quadratic scaling of the free energy holds is /-dependent and grows 
with increasing /. Indeed, the functionality acts as a prefactor that controls the strength 
of the logarithmic divergence of the potential at the origin. Formally, the mean-field 
approximation also becomes better with growing spatial dimension d, as the geometrical 
prefactor r"^"^ multiplying c(r) and — /5f (r) suppresses the small-r discrepancies of the 
two more efficiently. We call systems for which the mean-field idea holds only for the 
thermodynamics weak mean-field fluids. If one naively applies the strong mean-field 
relation, Eq. (pT]), to wieaA; mean-field fluids, one obtains results for the structure factor 
S{k) = [1 — pc{k)]^^ that can be seriously in error for finite fc-values. Only at = and 
at sufficiently high densities is it a reasonable approximation to set 5(0) = [l+pj3v{0)]~^. 
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Figure 2. Comparison of the mean- field approximation (dashed line) with the exact 
result (solid line) concerning the excess free energy density of the / = 32 star fluid. The 
slope of the straight line is 2, indicating the quadratic dependence of the excess free 
energy density on particle density. The arrow indicates the location of the crossover 



density px , above which the scaling of Eq. (15) holds with a relative error of less than 
10%. 



4. Clustering and reentrant melting 

In this section, we turn our attention to tlie pliase beliaviour of strong mean-field fluids. 
Two representatives of this class whose phase behaviour has been studied in detail are 
the Gaussian core model (GCM) of Eq. (Q) and the 'penetrable sphere model' (PSM) 
characterised by the interaction potential Vpsm{r) = eQ{a — r), with the Heaviside step 
function 6(x). Clearly, the PSM reduces to the hard sphere model for t = 0. 

The GCM has been the subject of extensive investigations by Stillinger et al. in the 



late 1970s |^ |3^, ^ The t = phase diagram of the model was calculated, 

showing the existence of two stable crystal structures, fee for low densities and bcc 
for high ones. In addition, a host of mathematical relations for the GCM has been 
established and on the basis of free energy estimates is has been postulated that the 
system displays reentrant melting behaviour at low temperatures. On the basis of 
simulation studies at selected thermodynamic points, a rough phase diagram of the 
GCM has been drawn A detailed study of the structural and phase behaviour of 



the GCM was carried out recently by Lang et al. There, it was indeed shown that 
for temperatures t > 0.01 the system remains fluid at all densities, whereas for t < 0.01 
reentrant melting is observed: increasing the density, the system first undergoes a fluid 
fee transition, followed by a structural fee — > bcc transition and at higher densities 
the bcc solid remelts, i.e., a bcc — > fluid transition takes place. The width of the solid- 
phase region grows with decreasing temperature. A structural signature of this unusual 
phase diagram in the fluid phase above t = 0.01 is an anomaly in the behaviour of the 
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liquid structure factor S{k). The height of its main peak first grows with increasing 
density and after achieving a maximum, it decreases again, reflecting the stability of 
the fluid beyond the reentrant melting. The Hansen- Verlet freezing criterion [53, 5^ 



was shown to be satisfied at both the freezing and the reentrant melting lines ||27|| . 
Moreover, it was found that at high densities not only the mean-field relation, Eq. (p^) 
is satisfied to excellent accuracy but also that the hypernetted chain closure (HNC) 



becomes quasi-exact |113|, The system becomes 'quasi-ideal' at those densities, 

meaning that the radial distribution function has the limiting behaviour g{r) — » 1. 
This, in conjunction with the mean-field property c(r) = —(3v{r) and the exact relation 
g{r) = exp[— /5f (r) — c(r) + g{r) — 1 + B{r)], forces the bridge function to obey the limit 
B{r) — > 0, hence rendering the HNC exact. 

The PSM was first studied in detail by by means of integral equation theories. 



computer simulations and cell-model calculations at small temperatures, t < 0.3 ||45 
In contrast to the GCM, no reentrant melting was found. Instead, the freezing line of the 
system appeared to persist at all temperatures, and cascades of clustering transitions 
in the solid were found, in which solids with multiple site occupancies are stable with 
increasing temperature and density. These findings were independently confirmed in 



a density-functional study of the low-temperature phase behaviour of the PSM |46 



Sophisticated integral-equation approaches at arbitrarily high temperatures revealed 
a loss of the solution along the 'diagonal' t = rj of the phase diagram again 
a feature pointing to an instability of the liquid by increasing density at arbitrarily 
high temperatures. Thus, the PSM and the GCM show completely different phase 
behaviours, although they are both bounded and non-attractive potentials. There is a 
cascade of clustering transitions for the former, enabling freezing at all temperatures, 
and a reentrant melting for the latter, associated with the inability to stabilise crystals 
above a certain critical temperature.|| 

The key in understanding these two very different types of behaviour lies in the 
strong mean-field character of these fluids and the associated expression for the fluid 
structure factor, Eq. (IT^) . If the Fourier transform of the pair interaction (j){k) has 
oscillatory behaviour (i.e., if 0(/c) becomes negative for some fc- values), then at the 
wavenumber k^^a where (j){ka) attains its most negative value, — |0(/c*(t)|, the liquid 
structure factor S{k^) will display a maximum. For any given temperature t, there 
exists then a density ps(t) such that ps|0(/^*o")| = t, causing a divergence of the fluid 
structure factor at k^, and marking a 'spinodal line' at this finite wavenumber. Thus, 
the fluid cannot be stable at all densities. As a matter of fact, freezing will take place 
before the spinodal line is reached. The PSM clearly belongs to this category, since the 
abrupt jump of the pair interaction fpsm('^) at r = cr causes long-range oscillations of 
the potential in Fourier space. By employing the Hansen- Verlet criterion, S{k^:) = 3, 

§ The use of the terminology 'critical temperature' here refers simply to the fact that above the said 
temperature freezing is impossible and should not be confused with its standard meaning in the realm 
of critical phenomena. There are no diverging thermodynamic quantities here, no universality and all 
free energies remain analytical functions in the neighborhood of the critical temperature. 
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it has been found |^T| that in the PSM freezing takes place at the 'diagonal' t = 77 on 
which the integral equation approach of Fernaud et al. breaks down. If, on the other 
hand, the Fourier transform of the potential, is a positive-definite, monotonically 
decreasing function of fc, Eq. (0) assures that at sufficiently high temperatures, where 
the mean-field approximation is valid at all densities S{k) is a monotonic function 
of k approaching rapidly the value S{k) = 1 with increasing k and being deprived of any 
peaks. The lack of peaks in the structure factor implies the lack of any tendency within 
the liquid towards spontaneous formation of spatially modulated patterns. Thus, the 
fiuid remains stable for all densities at sufficiently high temperatures. This, combined 
with the observation that at low temperatures and densities bounded potentials all 
become hard-sphere like and hence they must cause a freezing transition there, leads 
to a reentrant-melting scenario for such systems. Clearly, the GCM belongs to this 
category. Representative results for a particular family of strong mean-field fiuids and 
schematic phase diagrams can be found in Ref. pi] . 



5. Exotic crystal phases 

We now focus on weak-mean field fiuids, for which no simple criterion for their freezing 
behaviour can be established, since Eqs. (0) and ( ]T2| ) are not satisfied any more. The 
star-polymer fiuid characterised by the pair potential of Eq. (|) is a case in point. [j] The 
physical system of star polymers provides an excellent testbed for the investigation of 
the thermodynamic stability of more complicated crystals than usual fee- and bcc-lattice 
arrangements. 

The fcc-lattice is the one favoured by hard interactions, since it has the property 
of maximising the available volume and hence the entropy of the particles for a given 
particle density [Q. On the other hand, the presence of 'soft tails' in the potential has 
the effect of favouring the more open bcc-lattice, as was convincingly demonstrated for 
the case of the screened Coulomb potential (Yukawa interaction) arising in charge- 
stabilised colloidal suspensions |5^, ^ |5^. These two common lattices were 
considered for a long time to be the only 'candidates' in a search for stable crystals 
for given interatomic potentials. However, in modern colloidal science, new possibilities 
open up. It is technically possible to manufacture micelle-like particles featuring a 
hard core and a soft, fiuffy corona of grafted or adsorbed polymer chains, with the 
thickness L of the latter being much larger than the radius Rc of the former. Star 
polymers correspond to the case Rc ^ L; the theoretical arguments leading to the 
effective potential of Eq. (^) are in fact based on the assumption i?c — > 0. Under these 
physical circumstances, the effective interaction between the micelles is dominated by 
the ultrasoft repulsion between the overlapping, flexible coronas and not by the excluded 

II We consider here the case / > 10, for which indeed the interaction of Eq. (Q) holds, and not the low- 
functionality case for which the interaction of Eq. (^) is valid. The reason is that at low functionalities 
the stars do not freeze at any density and hence they are not an appropriate system for considering 
thermodynamically stable crystals. 
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volume interactions of the hard cores. Thereby, the requirement of maximising of the 
volume available to each particle does not play the decisive role any more. 

These properties manifest themselves in the phase diagram of star-polymer 
solutions. Due to the irrelevance of the temperature for this entropic interactions, 
the phase diagram was drawn in the (/, ?7)-plane by Watzlawek et al. [^, where 
rj = 7ipa^/6 and p = N/V is the number density of stars in the volume V. The phase 
diagram is shown in Fig. The fluid phase remains stable at all concentrations for 
f < fc = 34, a result that confirms and makes precise early scaling-argument predictions 
of Witten et al. For f > fc and at packing fractions 0.15 ^ ri < 0.70, the usual fcc- 
and bcc-crystals are seen to be stable, the former for larger and the latter for smaller 
functionalities. This is consistent with the fact that the effective potential of Eq. @ 
has a Yukawa decay length scaling as 1/ v7, hence large / is analogous to the strongly 
screened charge-stabilised colloidal suspensions. However, for rj > 0.70, unusual crystal 
structures appear. First, in the domain 0.70 < r/ < 1.10, a body-center-orthogonal 
(bco) crystal is thermodynamically stable. The bco-lattice is characterised by a body- 
centered, orthogonal conventional unit cell with three unequal sides and reduces to the 
bcc-lattice for ratios 1:1:1 between the sides and to the fee for ratios 1:1/ ■ 1 / 
|56| . The b co-lattices appearing in this region of the phase diagram feature strongly 



anisotropic unit cells with typical size ratios 1 : 0.6 : 0.3. Thus, these are crystals with 
coordination number 2. For packing fractions 1.10 ^f]^ 1.50, the diamond lattice with 
coordination number 4 turns out to be stable. Thus, we see that very open structures 
with their characteristically low coordination numbers are stabilised by the ultrasoft 
star-star potential. This feature has been attributed to the very slow divergence of 
the interaction as r — > 0, combined with its crossover to a Yukawa- form for r > a 
P, Bl, EM. Indeed, in such circumstances it may be energetically preferable for the 



system to have a small number of nearest neighbours at a small distance from any given 
lattice point than a large number of neighbours at a greater distance, as is the case for 
the optimally-packed fcc-lattice. 

In the fluid there exist clear structural signatures both for the topology of the 
phase diagram of Fig. ^ and for the variety of the crystal phases featured there. The 
reentrant transition from a fluid to a bcc-lattice and then again to a fluid, occurring 
for 34 < / < 48, is manifested in fluid structure factors S{q) that show a main peak 
that first grows with increasing density and then drops again as in the case of the 
Gaussian core model mentioned above. Moreover, once more the Hansen- Verlet freezing 
criterion [53, J4| was found to be satisfied on both sides of the freezing and reentrant 
melting line. The radial distribution function g{r) of the fluid at various densities, on 
the other hand, carries the signature of a local coordination that resembles that of the 
thermodynamically neighbouring solids. To demonstrate this, we show in Fig. ^ the 
function g[r) for star polymer fluids at / = 32, which are thermodynamically stable, 
at packing fractions rj = 0.80 and rj = 1.20. Comparison with Fig. ^ shows that the 
former corresponds to a state at the vicinity of the bco-phase and the latter to one at 
the vicinity of the diamond phase. Consider now the average coordination number z in 




Figure 3. The phase diagram of star polymers. The symbols denote simulation results 
for the pairs of coexisting densities and the lines are guides to the eye. Notice that the 
density gaps at the phase boundaries are very narrow. (Redrawn from Ref. [^.) 



the fluid phase, defined as 

z = Anp / r'^gir) dr, (1^ 



where r^am is the position for which g{r) has its first minimum and is indicated by the 
arrows in Fig. ^. For the two packing fractions shown we obtain z = 1.95 aX r] = 0.80 
and z = 4.03 at rj = 1.20. The first is very close to the coordination number Zbco = 2 
of the neighboring bco-lattice and the latter to -Zdiam = 4 of the diamond lattice. The 
fluid distribution functions contain local correlations that point to the ordering of the 
incipient crystal phases. 

The 5'(r)'s of the fluid above the overlap density, t] > 1.0, show in addition 
anomalous behaviour featuring two distinct length scales, as is clear from Fig. Kb). As 



analysed in detail in Ref. ||36|, two characteristics of the interaction are responsible for 
this behaviour: on the one hand, the existence of the crossover of the interaction of Eq. 
(1) at r = 0" from a logarithmic to an exponentially decaying form. And on the other, the 
ultrasoftness of the logarithmic potential, allowing the existence of fluids at arbitrarily 
large densities (for / < /c), a feature unknown for the usual interactions appearing in 
liquid-state physics and which are all 'perturbations' of the Hard-Sphere potential (e.g., 
Lennard- Jones, inverse-powers etc.) Thus, ultrasoft potentials carry unique structural 
signatures which should in principle be visible in scattering experiments from soft, 
polymeric fluids. 

A great deal of insight into the general physical mechanisms driving the stability 
of open structures in soft systems was gained through the recent work of Ziherl and 
Kamien [pTj They considered in full generality systems with particles composed of 
a hard core and long, deformable coronas and argued as follows. At any given density 




Figure 4. The radial distribution function of a star-polymer fluid of functionality 
/ = 32 at two different packing fractions, 77 — 0.80 [(a)] and 77 ~ 1.20 [(b)]. The arrows 
indicate the positions rmin that define the borderline of the first coordination shell in 
the fluid phase. 



above the overlap concentration, the coronas are forced to overlap and compress, which 
gives rise to an entropic free energy cost. The volume available to the coronas is fixed 
and equal to the difference of the total volume minus that occupied by the hard cores. 
Denoting by d the thickness of the coronal layer and A the total area of an imaginary 
membrane separating the compressed coronas, it then turns out that the product Ad 
is constant. As the free energy cost for the compression of the chains increases with 
decreasing thickness d^ it turns out that favourable phases are those for which the 
interfacial area A is minimal. Thereby, the problem reduces to that of determining the 
ordered arrangement of point particles that generates Wigner-Seitz (WS) cells having 
the smallest possible area for a given density. It is then conceivable that the fcc-lattice 
will be unfavoured, since its WS cell has a larger area than that of the bcc, for instance. 
In this way, Ziherl and Kamien established a beautiful connection of this problem with 
Lord Kelvin's celebrated question of determining the area-minimising partition of space 
for an arrangement of soap bubbles of equal volume |BU . 



Following these arguments, it then turns out that there exists yet another candidate 
phase that has an area even smaller than the bcc-lattice namely the Al5-lattice 
[l6l| . Self-assembled micelles of dendritic molecules with a particular architecture have 
been experimentally seen to crystallise into this phase |^ . The conventional unit cell of 
the Al5-lattice is shown in Fig. |^. It can be thought of as the cell of a bcc-lattice (dark 
points) decorated with 'dimers' (light points) running along the middle of the faces. 
The dimers are oriented parallel on opposite phases, forming thus columns through 
space. The orientation of the dimers lying on intersecting faces is perpendicular to one 
another, so that one third of all dimers lies along each of the three Cartesian directions 
in space. The dimer length is a/2, where a is the edge length of the cube, and it is placed 
symmetrically along the face, i.e., the distance of any monomer to its nearest edge is 
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a/4. The A15- lattice is not a Bravais lattice; it can be constructed as a simple cubic 
(sc) arrangement with an eight-member basis, thus it contains 8 sites per conventional 
cell. Its WS-cell is a Goldberg decatetrahedronpj] consisting of two hexagonal and 12 
pentagonal faces fS^ . 

Semi-quantitative calculations on the stability of the Al5-lattice were carried out 
by Ziherl and Kamien ||5^ , using a model 'square-shoulder' potential within a simplified 
cell model. Narrow regions in thermodynamic phase space were found, in which the A15- 
lattice was stable but this finding is uncertain in view of the approximations involved 
and the limited extent of this region. To investigate this question in more detail, we 
have employed extensive lattice-sum calculations for the star-polymer system, using the 
pair potential of Eq. and extending both the set of candidate lattices and the region 
of densities we looked at. We compared between the sc, diamond, bco (which includes 
the fee- and bcc-lattices as special cases) and Al5-lattices in the regions < ?7 < 2.50 
and 32 < / < 256, at selected arm numbers / to be shown below. The lattice sums 
were performed by fixing the particles at the prescribed lattice positions and keeping 
them frozen there, i.e., no thermal fluctuations (harmonic corrections) were taken into 
account. This approach reproduces very well the solid part of the phase diagram of Fig. 
^: the free energy of the star-polymer crystals turns out to be dominated by the lattice- 
sum term, the corrections to it from the particle oscillations as well as the entropic 

% We prefer the term decatetrahedron for a polyhedron with 14 faces instead of the term 
tetrakaidecahedron, often employed in the literature. The latter, inspired from ancient Greek, literally 
means 'four-and-ten-faced polyhedron', whereas the former, consistent with modern Greek, has the 
much more logical translation 'fourteen- faced polyhedron'. One of us (CNL) believes that modern 
Greek language deserves a fair chance against its classical predecessor, which has influenced scientific 
terminology for quite some time. 
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contribution from the same playing only a minor role. In this way, we are of course 
unable to compare the solid free energies with those of the fluid, therefore no prediction 
about melting can be made. However, for large enough /, the interaction is steep enough, 
so that the system will be deflnitely be in a crystalline phase for which the lattice sums 
provide an accurate prediction of the most stable structure among the candidates. 

For the bco-phases we minimised the lattice sums with respect to the two size 
ratios, ri = b/a and r2 = c/a between the edge lengths, a, b and c at any given density. 
Without loss of generality, we assume in what follows that a is the longest of the three 
edges, thus < ri,r2 < 1. The minimised bco-energy was then compared with the 
energies of all other lattices; the one with the smallest lattice energy per particle wins. 




I I 

0.0 0.5 1.0 1.5 2.0 2.5 



Figure 6. The zero-temperature phase diagram of star polymer solutions, obtained 
after the minimisation of lattice sums for particles interacting by means of the potential 
of Eq. {4). The phase denoted bcoi in this figure is the same as the one denoted bco 
in Fig. B but now it has to be distinguished from the additional bco-phases showing 
up at higher densities and named bco2, bcos and bco4. 



The 'zero-temperature' phase diagram^ obtained this way is shown in Fig. ^. First, 
we note that the phases being stable up to = 1.50 are precisely those also seen in the 
flnite-temperature phase diagram of Fig. |] and also that the phase boundaries based on 
the lattice sums agree very well with the ones at finite temperatures. The Al5-phase 
does not alter the hitherto explored part of the phase diagram. However, for r] > 1.50, 
a host of new phases and of transitions between those show up. The Al5-phase is stable 
at the high-density part of the phase diagram, confirming thus explicitly the prediction 

'Zero-temperature' is just a convenient way to refer to the assumption of frozen particles at the 
lattice sites and does not refer to the real temperature T of the system. The latter is an irrelevant 
thermodynamic variable because the entropic interaction of Eq. is proportional to ksT and the 
thermal energy is the only energy scale of the problem. 
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of Refs. 1^ and |0 that this ph suitable candidate at the high concentrations 

of ultrasoft particles. Nested between the stability domain of the Al5-lattice and the 
diamond lattice, four new bco-phases and (iso) structural transitions between those show 
up, having the following characteristics. 

The phase denoted bco2 has size ratios ri = 1 and r2 = 0.55, the former being 
constant at all densities and functionalities and the latter showing very weak variation, 
see also Fig. 0. Hence, the unit cell of the bco2-phase is anisotropic only in one Cartesian 
direction and has a wide basis and height that is smaller than the base edge length. 
Accordingly, the coordination number of this phase is 2. The bcos-phase, dominating 
at high functionalities, has size ratios ri = r2 = 0.3 which also show very little variation 
with rj and /, see Fig. ^ Similarly to the bco2-phase, therefore, the cell of the bcoa-phase 
also has anisotropy in only one Cartesian direction. Contrary to it, however, the height 
of the cell is now larger than the edge length of the base and therefore the number of 
nearest neighbours is 4. The phase transition bco2 bcoa is first order: as can be seen 
in Figs. |^(c) and (d), the size ratios jump abruptly from the values (1, 0.55) of bco2 to 
the values (0.3, 0.3) of bcoa. These are the only stable bco-phases in this part of the 
phase diagram as long as / > 56. 

For functionalities 32 < / < 56, two more bco-phases show up. First, there is 
a narrow region occupied by the bco4-phase, which has all three edge lengths of its 
unit cell different and is thus similar to the bcoi-phase discovered before. As seen in 
Figs. 13(a), (b) and (c), these ratios evolve from the values ri = 1 and r2 = 0.55 of the 
bco2-phase towards the values ri = l/\/2 and r2 = l/\/2 of the fcc-phase. 

The order of the transitions between those phases has been investigated numerically. 
Within the limits of accuracy of our numerical code, we have found that for / = 32 and 
/ = 40 the transition bco4 — > fee is first order, i.e., the size ratios of the bco4-phase 
jump with density to the ratios ri = r2 = 1/ of the fcc-phase abruptly. This is seen 
in Figs. |^(a) and Fig. 0(b). For / = 48, the transition is second-order with the bco4-size 
ratios evolving to the fcc-ones smoothly, see Fig. |^(c). The nature of the transition bco2 
— > bco4 is also mixed. For / = 32, [Fig. |^(a)], a finite jump of the values of the ratios 
was found, although the step size in changing the packing fraction was made as small as 
5 X 10~^. Thus, we characterise this transition at / = 32 as first-order. For / = 40, [Fig. 
0(b)] this transition appears to be second-order. Hence, we conclude that there must 
be a line of second-order transitions terminating at a tricritical point between / = 40 
and / = 32, (for the bco2 bco4 transition) and similarly a tricritical point between 
/ = 48 and / = 40 (for the bco4 — > fcc-transition) to be succeeded by lines of first-order 
transitions. 

All four bco-phases and the three transitions among them, bco2 bcos hco^ 
fee, appear only in a very narrow /-range around / = 48, see Fig. 0(c). The line of 
second-order transitions bco2 — bco4 meets the lines of first-order transitions bco2 — *■ 
bco3 and bcoa —>■ hco^ at a critical endpoint located at / slightly less than 48. Similarly, 
the line of second-order transitions bco4 fee also terminates at a critical endpoint, 
meeting the first-order lines bcos ^ A15 and fee — A15. It is an intriguing phenomenon 
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Figure 7. The optimal size ratios of the various bco-phases occurring in the phase 
diagram of Fig. ^ between the diamond and the A15-phase, within their domains 
of stabiHty. Also shown is the characterisation of those phases. The four different 
functionalities are indicated on the plots. 



that such a richness in the stable crystal structures and in the nature of the transitions 
among them occurs as the result of a simple, spherically symmetric interaction and this 
points to the many surprises of ultrasoft potentials and their tendency to produce open, 
exotic structures. At the same time, it must be pointed out that we do not expect the 
solid phases occurring for / < 48 to survive the competition with a fluid. Indeed, as 
can be seen from Fig. ^, the critical functionality fc below which no solids are stable 
increases with increasing density. The bcc-phase is extinguished by the fluid for / < 34, 
whereas the bco- and diamond phases for / < 44. It is therefore anticipated that the 
bco2-, bco4- and fcc-phases seen in Fig. ^ for / < 48 will be wiped out by the fluid 
there. However, by arbitrarily increasing / one can always reach a domain where the 
fluid will be beaten by the crystal and hence the bco2 —>■ bcoa —>■ Al5-transitions will 
be there also after a finite-temperature calculation. In particular, the Al5-phase, on 
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whose stability it has been speculated for some time, has now been proved to be indeed 
the most stable phase at sufficiently high densities among the candidates considered. 

The logarithmic- Yukawa potential employed in our study has been shown to 
describe well the effective interactions between star polymers in a good solvent for a 
wide range of concentrations. However, at the region of stability of the Al5-crystal, 
the pair potential description is not expected to be particularly accurate. Many-body 
contributions are expected to become important there |6^. Moreover, the Yukawa 
tail of the potential, describing the interactions of the outermost Daoud-Cotton blobs, 
should be absent since the compressed coronas there are deprived of the outermost blob 
structure of the isolated stars. Thus, in this respect, the logarithmic- Yukawa potential 
has to be looked upon rather as a toy model. Nevertheless, the physical characteristic 
driving the transitions discovered above is the ultraslow divergence of the logarithm 
which, in the neighborhood of the average particle distance can be locally expanded as 
a ramp-like potential and the finer details of the interaction should become irrelevant. 



6. Summary and conclusions 

We have shown that "ultrasoft interactions" arise naturally from coarse-graining 
procedures for a broad range of soft matter systems. Besides greatly simplifying the 
statistical mechanics of these complex systems - once the interactions are derived, all 
the well known tools of liquid state theory can be applied to calculate correlations and 
phase behaviour - they also lead to new phenomenology. Signatures of these ultra 
soft interactions include anomalous fluid correlations, reentrant melting as well as the 
stabilisation of exotic, open crystal structures. In contrast to their atomic counterparts, 
soft matter systems can therefore stabilise such crystals without the presence of angle- 
dependent, anisotropic potentials: radially symmetric, ultrasoft interactions are quite 
sufficient. Thus, a new mean fi.eld fluid paiadigm can be established, which goes beyond 
the usual prototype for classical fluids, the hard-sphere model. The latter, being always 
dominated by packing effects, tends to favour close packed structures. Exotic crystals 
with unusual ordering have been observed in hard-sphere-like suspensions ||6^, ^ but 
that case refers to binary mixtures whose phase behaviour is indeed much richer than 
that of their one-component counterparts. 
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